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Abstract 

Living communities can be considered as complex systems, thus a fertile ground for studies 
related to their statistics and dynamics. In this study we revisit the case of the rhythmic applause 
by utilizing the model proposed by Vazquez et al. [A. Vazquez et al., Phys. Rev. E 73, 036127 
(2006)] augmented with two contradicted driving forces, namely: Individuality and Companionship. 
To that extend, after performing computer simulations with a large number of oscillators we propose 
an explanation on the following open questions (a) why synchronization occurs suddenly, and b) 
why synchronization is observed when the clapping period (Tc) is 1.5 • < Tc < 2.0 • (T^ is the 
mean self period of the spectators) and is lost after a time. Moreover, based on the model, a weak 
preferential attachment principle is proposed which can produce complex networks obeying power 
law in the distribution of number edges per node with exponent greater than 3. 
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I. INTRODUCTION 



Periodic phenomena are of high abundance in nature, what makes them attractive to 
physicists is the fact that they can serve as prototype complex networks, thus the impor- 
tance on understanding their statistics and dynamics. To give an idea about the universahty 
of the phenomenon we mention, among other systems, the: Pteroptyx malaccae fireflies [1], 
secretory cells [2], synchronous firing neurons [3j, and rhythmic applause [4j. In particular, 
the study of the mechanism that leads to self organization of biological systems through syn- 
chronization is of major importance, since it provides enough information as to understand 
the dynamics of the interactions among the members of the system. 

To address the problem of synchronization we used the coupled oscillators approach. It is 
believed that Huygens ideas about the synchronous move of wall hanged pendulums initiated 
the studies towards the explanation of this phenomenon. It is noted that coupling between 
identical oscillators is trivial since it can be achieved by a phase-minimization procedure. 
This is not the case in biological systems since the oscillators can be considered anything 
but equivalent. In such diverse systems depending on the strength and type of interaction 
as well as the dispersion of the self periods, synchronization might or might not appear. To 
illustrate the problem we refer to the synchronous clapping of an audience after a spectacle. 
The audience is comprised by a large number of spectators each one characterized by her/his 
temper, enthusiasm, acoustic behavior, thus can be considered to applause in her/his unique 
self period. 

A first approach to the rhythmic applause indicates that the integrate-and-fire-type model 
[HI E] is the most promising to understand the synchronization mechanism. However, this 
model doesn't take into account memory effects which are crucial for the synchronization 
to occur, thus the continuous phase coupling in the Kuramoto model [7| has been used 
to understand the underlying dynamics [H]. On the other hand, a recent study brought 
fourth a strong evidence that human dynamics of many social, technological and economic 
phenomena can be modeled by a decision based queuing process |9] in which individuals 
execute a task from a list of pending tasks based on a given priority. 

The purpose of the present paper is to show that application of the aforementioned de- 
cision based queuing process is capable of modeling the dynamics of the rhythmic applause 
and especially: i) why synchronization occurs suddenly, ii) why synchronization is achieved 
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in a period almost two times the average self period of the spectators, and iii) why synchro- 
nization appears and disappears several times during the applause. 



II. RHYTHMIC APPLAUSE AND THE KURAMOTO MODEL 

The Kuramoto model has been used to model the physics of the rhythmic applause. In 
this model we deal with N-oscillators each of them described by its (f)j phase. The oscillators 
have a g{uj) distribution of their ujj self frequency (where g{uj) is considered to be a normal 
distribution). Every rotator interacts with all other ones via phase-difference-minimizing 
terms via the following formula 

N 

N 



i=l 

The N coupled differential equations describing the over-damped oscillator dynamics are 

i=l 

Mathematically the synchronization level will be characterized by an order parameter, q, 
defined at any time moment as 

1 ^ 
N ^ 

The maximal possible value q=l corresponds to total synchronization, the case < g < 1 
to partial synchronization, while for q=0 there is no synchronization at all in the system. 
In the N ^ oo thermodynamic limit of the equilibrium dynamics ^ oo, so initial 

transient effects are lost) Kuramoto and Nishikawa proved the existence of a K^. critical 
coupling. For a Gaussian distribution of the oscillators' natural frequencies, characterized 
by a dispersion D, they got 
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(3) 



K.= \I^,D (4) 



For K < Kc, the only possible solution gives q=0 (no synchronization) while for K > 
a stable solution with g 7^ appears. Thus, the main result is that for a population of 
globally coupled nonidentical oscillators a partial synchronization of the phases is possible 
whenever the interaction among oscillators exceeds a critical value. 

There are for four major drawbacks of this model when applied to rhythmic applause: 



1. The uniform one to one coupling between all the oscillators. The spatial distribution 
of spectators prohibits this assumption. 

2. The coupling constant should be K > Kc. Since i) synchronization is achieved at a 
period ^ 2 ■ (with T, the mean self period of the spectators), and ii) the doubling 
of the period leads to the half-doubling of the dispersion coefficient (D), the following 
inequality 



should hold. We note that D is strongly related to behavioral characteristics of the 
spectator (temperament, enthusiasm), thus unrelated to K which, in turns, is related 
to factors such as the distribution of the sound in the room, reflections of the sound etc. 
On the other hand the "generalized relaxation oscillators" approach is characterized, 
again, by the need of a global coupling among the oscillators that can not be guaranteed 
in networks of living organisms. 

3. Usually at the beginning of the applause there is a long "waiting" time without any 
synchronization, and with no increase in the order parameter. Partial synchronization 
evolves suddenly after that and achieves its maximal value in a short time. This 
should not be the case if we are in the K > Kc limit. One would expect in this limit a 
continuous increase in the order parameter right from the beginning of the applause. 

4. why is synchronization already achieved lost after a time, and why might it reappear 
again? Loss of synchronization should not happen in the K > Kc limit. 

III. THE DECISION BASED QUEUING PROCESS 

Recently, Vazquez et al. P proposed a decision based queuing process (DBQP) able 
to reproduce the distribution of the possible action of humans in particular cases. In this 
model, all the possible actions are forming a list and the acted one is chosen by applying 
criteria (e.g., priority) related to the specific individual, that is the emotional factors are 
playing an equally important role. 

Applying the above DBQP model in the case of the rhythmic applause, the spectator 
she/he has to chose among the following actions (tasks): 




(5) 
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1. to clap at the frequency that fulfill her/his satisfaction (this an expression of Individ- 
ualism), or 

2. to try to change her/his pace as to follow the clapping frequency of the synchronized 
spectators (this an expression of Companionship). It is obvious that this holds true 
only in the presence of such a group of people. 

Every time that a spectator claps her/his hands she/he makes a decision between (1) 
and (2), depending on the priority given to each of them. The priority is immediately 
related to the reached level of satisfaction reflected to either the need of expressing her/his 
Individualism or her/his Companionship. We name the flrst as Iprio- The fluctuation of this 
parameter with the time is ruled by the level of its saturation: when it is saturated then Iprio 
reduces with time, whereas increases in the other case. Two boundaries control this behavior, 
namely the High Saturation Threshold (HST) and the Low Saturation Threshold (LST). 
When Iprio climbs over HST, then the spectator's need to express her/his Individualism 
is considered to saturate, thus the decay of Iprio- On the other hand when Iprio fall under 
LST, then the need to express the Individualism comes again in the foreground, thus the 
Iprio starts to increase. This mechanism is depicted in Fig. [l] 

A model system of a square grid with n x n = n"^ spectators is chosen and it is assumed 
that the sound produced in position, {i,j) is heard weakened, proportional to the distance 
from the source, in neighboring regions. Thus, the intensity of the sound in position {i,j) 
at time t is: 

n n 

si^.j. ^) = E E ^) ■ - ■ - ^■') (6) 
i'=i j'=i 

where K is the sound weakening factor (a 5-point derivative kernel [ID]), and 6{i',j',t') 
takes values 1 and if at time t the spectator at position {i',j') claps or not her/his hands, 
respectively. We assume that at time t the sound intensity in the {i,j) position is S{i,j,t) 
and S{i,j, t) = 1, meaning at that time the spectator at position {i,j) claps her/his hands. It 
can be shown that the S{i, j, t) could give us information on the total number of synchronized 
clapping spectators. At time t' when 6{i,j,t') = but S{i,j,t') > S{i,j,t) the sound heard 
at position {i,j) is stronger to that at time t when (i, j)-spectator also acted. This emerges 
the need of expressing her/his Companionship. The greater the S{i, j, t') / S{i, j, t) the higher 
the priority given to this action, namely Cprio- 
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At any given moment t a random number 5" is generated: if S < {Iprio/ilprio + C'prio)) 
then the action (1) is chosen else the action (2). 

Let us see what happens when the action (2) is chosen. If t is the present time and the 
time (before t) when S{i,j,ta) = 1 and tb the time (after t) when 6{i,j,tb) = 1 again, then 
h = 'ta + Ti,j, with Tij the self-period of the spectator (z, j) (see Fig. |2]). 

The following three cases might occur: 

• if t is closer to ta, the spectator increases her/his self-period hoping that she/he will 
be synchronized with the next clapping. 

• if t is closer to tb, the spectator lowers her/his self-period hoping that she/he will be 
synchronized with the next clapping. 

• if t is closer to {ta + tb)/2, the clapping at that time will result to an increment of the 
background noise that is strongly related to the synchronization. To that extend, the 

wishes not to increase the background noise thus suspends her/his clapping at 
time tb (we note that this could happen only once). 

IV. EXPERIMENTAL RESULTS - DISCUSSION 

A statistical ensemble resembling a grid of 625 (25 x 25) spectators was simulated. We 
begin the simulation with self-periods following normal distribution with average value 50 
and a = 10. We assume that immediately after the clapping starts, spectators feel a great 
need for expressing her/his enthusiasm or Individualism, thus Iprior is high. For this reason 
spectators retain their self-periods of clapping and consequently there are no signs of syn- 
chronization. When time passes Iprior gets lower resulting to the forming of random, at the 
beginning, clusters, that give rise to the expression of Companionship. Immediately after 
the forming of the first cluster an abrupt increase of the average period is observed as it 
shown in Fig. |3} This is in agreement with results presented in [U [S] . 

The critical point tc of this abrupt change of period {tc ~ 2000) is characterized by a 
domination of the Cprio over the Iprio and for this to happen, Iprio has to be low enough, 
thus the incoherent clapping before synchronization occurs. After some time, Iprio starts to 
increase because it reaches the LST. Spectators are starting to express again themselves by 
clapping at their eigenfrequency resulting to the loss of synchronization. Later, Iprio starts to 
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decay because it reaches the HST and synchronization might occur again. This mechanism 
might be repeated several times during the clapping period as it is found experimentally [H] . 

Figure |4] shows the sound level at the position (12,12) of the audience. At the critical time 
tc, synchronization is visible as pronounced spikes in accordance to experimental findings 
[8] . The dominant period after time tc can be evaluated by counting the spikes resulting to 
a value of ~ 100 which is in agreement with the findings of Fig. [3j 

Moreover, it seems that the mean period of the spectators is a function of LST: the 
maximum mean period recorded in the simulation vs. the LST is depicted in Fig. [5j The 
error bars are proportional to the background sound, which is a measure of synchronization 
as it is explained in [1]. 

Finally, it is interesting to follow the dynamics of the abrupt change of the mean period 
shown in Fig. [3} Initially, lack of synchronization is due to the high value of the Iprio as it was 
explained before. After some time (the value of Iprio has been decreased), random clusters 
of synchronized spectators may be formed producing locally a high sound level. Spectators 
near those clusters either slowly start to take action to synchronize with them or select to 
suspend some of their claps. In the latter case the background noise is decreased (since some 
claps are missing) causing Cprio to increase suddenly. This effect acts as an avalanche making 
the other spectators to act faster to synchronize with the already formed clusters since their 
sound will now be more evident. Hence, it is expected that the aforementioned mechanism 
has a stronger effect on those spectators who clap at higher frequencies, since these spectators 
contribute more to background noise. For this reason, low clapping frequencies are sustained 
and the dominant period during synchronization is increased in comparison to the mean self 
period of the spectators. 

Our proposal was further tested for its validity against the Barabasi-Albert model [HI [12] 
that introduces the "preferential attachment" concept in the interpretation of the dynam- 
ics of large and complex networks. We formulate the model by setting up the following 
differential equation: 



where kj, ki is the degree of the respective node (the number of edges starting from node 
j or i), t is the time step, m is the number of links that Nth node will establish, — 1 
the number of the existing nodes of the cluster, and a is a parameter that takes values as 




(7) 
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< a < 1; when a = 0(a = 1) it is the C ompanionship{Individuality) that controls the 
direction of the added hnks. In other words the left term from the right hand part of Eq. [7] 
counts for the expression of Individuality as: "links will be placed randomly", whereas the 
right term counts for the expression of Companionship as: "links will be directed to the 
node with the highest degree" . The above Eq. [7] can be produced by the following queuing 
process; when a new node is added to the network, the m-edges that carries are placed 
one by one as follows: The first node is placed randomly with probability pi or following 
the "preferential attachment" mechanism with probability 1 — pi- Next, the second edge 
is placed following the same scheme but with probability p2, and so on. It can be shown 
that a = ^ Pi Eq. [7] holds. The solution of this differential equation leads to the 
following: 

where /3 = 1^ and ti the moment of insertion of the new node. Hence, 



l-(3 fty l-2f3 

-ml— >> — - — m (9j 



then the Eq. [8] is reduced to: 



(3 \t 

That is the distibution of degrees is given as: 



k,{t)=^—^m(U (10) 



P(t) ^ 2 ( (11) 



thus, the network is evolved by following a power law defined as 



7 = 1 + ^ or 7 = 1 + > 3, for a ^ 1 (12) 

fj 1 — a 



In the extreme when a = 1 the model reduces to random network [13], while for a = 
the model coincides with the Barabasi- Albert model [11]. The above described formalism 
enables us to produce undirected scale free networks with 3 < 7 < oo. An illustrative 
example comes from the preferential attachment in sexual networks. In this particular case 
the value of 7 = 3.4 is the largest observed [121 El [E] resulting to a = 16.7% or an 
expression of Companionship which in this case is translated to: "the target sexual partner 
is the one with the highest number of sexual links". This argument further supported by 
the fact that although preferential attachment in sexual networks has been experimentally 
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established ^6], the individual heterogeneity in the inclination to find new partners was 
found essential to model experimental data [T3] . 



V. CONCLUSIONS 

In summary, by augmenting the Vazquez et al. [S] model with two contradicted driving 
forces we were able to simulate the abrupt nature of synchronization. In addition we have 
shown that since in the first moments after the performance most spectators applause at 
their own pace {Iprio is large) there is no synchronization. Later, when the value of Iprio has 
been decreased, spectators may act to express their Companionship. Random clusters of 
partial synchronized spectators cause an abrupt increment of the Cprio of the neighboring 
spectators, which in turns leads to "global" synchronization. However, after some "synchro- 
nized" period the Jp^io could be enhanced again, resulting to lose of synchronization. Last 
but not least the expression of Individualism does not characterize communities of the East- 
ern European countries where Companionship is the main driving force, thus the easiness to 
achieve synchronization, whereas the opposite is observed in Western countries communities 

HI- 
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List of Figures 

Fig. 1 The priority Iprio of spectators to express their Individualism as a function of time. 
When Iprio chmbs over HST then the need is considered saturated and Ip^o starts to 
decay. The opposite behavior for Ip^o is observed when it reaches LST, and thus the 
need comes again to the foreground causing Iprio to increase. 

Fig. 2 Different moments during the decision making interval. 

Fig. 3 The average period of an audience of 625 spectators in a rectangular 25 x 25 grid. 
The self-periods follow a normal distribution with mean value 50 and cr = 10. The 
simulation parameters are: HST = 0.9, LST = 0.05, diprio/dt = ±0.001. 

Fig. 4 The sound level at position (12,12) of the audience of the same simulation as in Fig. 
[3] Synchronization after the critical time is indicated by the recorded spikes. 

Fig. 5 The maximum mean period recorded just after the critical time tc, of the same 
simulation as in Fig. |3j vs. the level of the Low Saturation Threshold LST. The error 
bars are proportional to the background noise which is considered to be a measure of 
lack of synchronization. 
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